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Abstract 

We have reevaluated the hadronic contribution to the anomalous magnetic moment of the muon 
(g — 2) and to the running of the QED fine structure constant a(s) at s = M|. We incorpo- 
rated new data from hadronic r decays, recently published by the ALEPH Collaboration. In 
addition, compared to previous analyses, we use more extensive e + e~ annihilation data sets. 
The integration over the total hadronic cross section is performed using experimental data up to 
40 GeV and results from perturbative QCD above 40 GeV. The improvement from r data con- 
cerns mainly the pion form factor, where the uncertainty in the corresponding integral could be 
reduced by more than a factor of two. We obtain for the lowest order hadronic vacuum polariza- 
tion graph a|f d = (695.0 ± 15.0) xl0~ 10 and Aa^Aff) = (280.9 ± 6.3) x 10" 4 using e+e~ data 
only. The corresponding results for combined e + e~ and r data are aj} ad = (701.1 ± 9.4) x 10 -10 
and Aa[^ d (M|) = (281.7 ± 6.2) x 10 -4 , where the latter is calculated using the contribution from 
the five lightest quarks. 
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Introduction 

The anomalous magnetic moment a M = (g — 2)/2 of the muon is experimentally and the- 
oretically known to very high accuracy. In addition, the contribution of heavier objects 
to a M relative to the anomalous moment of the electron scales as (m M /m e ) 2 ~ 4 x 10 4 . 
These properties allow an extremely sensitive test of the validity of QED and of addi- 
tional contributions from strong and electroweak interactions. The present value from the 
combined /x + and fi~ measurements [jlj, 

a M = (11 659 230 ±85) x 1(T 10 , (1) 

should to be improved to a precision of at least 4 x 1CT 10 by a forthcoming Brookhaven 
experiment (BNL-E821) @j. 

It is convenient to separate the prediction from the Standard Model (SM) a^ M into its 
different contributions 

af = af D + aj ad + a^ k , (2) 

where a^ ED = (11 658 470.6 ± 0.2) x 10~ 10 is the pure electromagnetic contribution (see @ 
and references therein), a^ ad is the contribution from hadronic vacuum polarization, and 
a" eak = (15.1 ± 0.4) x 10~ 10 |3], |4], || accounts for corrections due to the exchange of 
the weak interacting bosons up to two loops. Using the recent analysis of S. Eidelman 
and F. Jegerlehner || that found aj^ ad = (702.4 ± 15.3) x 10~ 10 , and applying fourth 
order corrections due to the exchange of additional photons and electron or quark loops 
(summarized by T. Kinoshita et al. to (—4.1 ± 0.7) x 10 -10 ) one finds 

a^ M = (11 659 184 ± 16) x 10~ 10 . (3) 

Comparing the errors of the respective contributions to a™ reveals that its total un- 
certainty is clearly dominated by the leading order vacuum polarization correction a^ ad , 
originating from a quark-loop insertion into the muon vertex correction diagram as shown 
in Figure [I]. 

In this letter, we present a new evaluation of the hadronic vacuum polarization con- 
tribution to a M and also to the running of the QED fine structure constant a(s) from low 
energy to the mass of the Z boson. In addition to using the complete and in comparison 
with previous analyses slightly enlarged experimental information on e + e~ annihilation 
data, we incorporate new data from hadronic r decays |]] which provide a more precise 
description of the hadronic contributions at energies less than 1.5 GeV. We bring attention 
to the straightforward and statistically well-defined averaging procedure and error prop- 
agation used in this paper, which takes into account full systematic correlations between 
the cross section measurements. We also stress our careful treatment of unmeasured final 
states which are bound via isospin constraints. 

1 Hadronic Vacuum Polarization in 7 and W propa- 
gators 

As QCD is a non-Abelian theory with massless gauge bosons, its perturbative expansion 
at low energies is not well-behaved and non-perturbative effects lead to currently un- 
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predictable long distance resonance phenomena in quark interactions. Fortunately, cross 
sections measured in e + e~ annihilation and spectral functions from r decays provide an 
experimental access to the hadronic vacuum polarization: from unitarity, the hadronic 
cross section of e + e~ annihilation is related to the absorptive part of the vacuum polar- 
ization correlator via the optical theorem. 

Similarly, hadronic spectral functions from r decays are directly related to the isovec- 
tor vacuum polarization currents when isospin invariance (CVC) and unitarity hold. For 
this purpose we have to worry whether the breakdown of CVC due to quark mass effects 
(Tn u 7^ ma generating d p J^ ~ (m u — rrid) for a charge- changing hadronic current J M 
between u and d quarks) or unknown isospin-violating electromagnetic decays has non- 
negligible contributions within the present accuracy. Recent estimates of CVC predictions 
of r branching ratios into vector final states |9j however show good agreement within 
about 5% experimental accuracy over the full range of exclusive vector hadronic final 
states. Expected deviations from CVC due to so-called second class currents as, e.g., 
the decay t~ — > 7i^f]u T where the corresponding e + e~ final state 7i°r] (C=+l) is strictly 
forbidden, have estimated branching fractions of order of (m u — rrid) 2 — 10~ 5 [III]], while 



the experimental upper limit amounts to B(r — > ir~r]is T ) < 1.4 x 10 |TT| with 95% 
CL. Another classical test is the pion /5-decay, yielding a sensitivity to CVC violation 
of 3%. The CVC hypothesis relates the isovector, vector matrix element of the decay 
7r~ — > 7r°e~z/ e to the electromagnetic form factor of the pion. No deviation between the 



CVC branching ratio and the experimental result has been found [12, 11 



An estimate of a possible CVC violation can be obtained in the nir final state which 
is dominated by the p(770) resonance. SU(2) symmetry breaking can occur in the p^-p 
masses and widths caused by electromagnetic interactions. Hadronic contributions to the 
p ± -p° width difference are expected to be much smaller since they are proportional to 
(m u — rrid) 2 . The different electromagnetic contributions to the width are listed in Table [TJ: 
radiative transitions introduce a negligible effect, while the dominant contribution comes 
from the n^-n mass difference. A recent theoretical analysis |TJ]] indicates that the p^ 
and p° masses are in fact equal within 0.5 MeV/c 2 : m p ± — m p o = —(0.15 ±0.55) MeV/c 2 . 
This prediction has been verified in the ALEPH analysis of the r vector spectral func- 
tions || with the result: m p ± — m p o = —(0.0 ± 1.0) MeV/c 2 . Since the n~ and it mass 
difference is known experimentally |ll| to be m^± — m^o = (4.5936 ± 0.0005) MeV/c 2 
and understood theoretically |15| to be almost completely from electromagnetic origin, it 
is expected that the total p ± and p° widths should be different even in the limit where 
hadronic interactions are S*?7(2)-invariant (this includes the chiral limit). An additional 
point concerns contributions from photon bremsstrahlung. While the infrared divergences 
in p —>■ 7T7T7 decays vanishes when including the vertex correction graphs, finite terms are 
expected to contribute differently to the widths of the charged and the neutral p. The 
corresponding bremsstrahlung graphs have been calculated in Ref . [13| . The width contri- 
butions from finite terms to both the p^ and the p° turn out to be negative. The estimate 
of the width difference given in Tab. [I] assumes finite contributions to the widths from 
loop corrections to be small. 

The total expected SU(2) violation in the p width is finally computed from the above 
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Table 1: Expected CVC violation from electromagnetic interactions in p ± -p° decays. 
considerations to be 

Tp± ~ Tp ° = (2.8 ± 3.9) x 10~ 3 , (4) 

where the error comes essentially from the estimate of the p ± -p° mass difference and of 
the 7T7T7 contribution. This effect introduces corrections for a^ ad and the running of a(s) 
when including r data (see Section^). 

From a more qualitative point of view, one should keep in mind that in this analysis 
the CVC hypothesis is applied at a very low energy scale where the absorptive parts of 
the matrix elements are largely dominated by non-perturbative QCD which are expected 
to factorize from their respective W or 7 excitation. 

However, electroweak radiative corrections must be taken into account. Their domi- 
nant contribution comes from the short distance correction to the effective four-fermion 
coupling r~ — > (du)~ v T yielding the fractional logarithmic term (a(M T ) /2tt) (4 ln(Mz/M T ) 
±5/6) [16| which vanishes in leptonic r decays. The short distance correction can be 
absorbed into an overall multiplicative electroweak correction Sew introduced in the 
definition of the spectral functions used in Ref. 0. The resummation of higher-order 
electroweak logarithms using the renormalization group yields [T7| , 18 ] 



while remaining perturbative electroweak corrections are of order a n (M T ) ln n (Mz/M T ) 
which is safe to ignore. The subleading non-logarithmic short distance correction, cal- 



culated to order O(a) at quark level: 5a(M T )/127r ~ 0.0010 [16] turns out to be small. 
Additional long-distance corrections are expected to be final state dependent. They have 
only been computed for the r~ — >tt~ v t decay leading to a total radiative correction of 



2.03% [B| which is dominated by the leading logarithm from the short distance contribu- 
tion. The evaluation of @ neglects radiative corrections from additional gluon exchange 
between the final state quarks. 
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had 

Figure 1: Leading order hadronic vacuum polarization contribution to a M . 



To be safe |2"0|] , the uncertainty of Sew in © is estimated to be A Sew — 0.0040, 
which is taken into account in the CVC cross section prediction from r decays. 

Muon Magnetic Anomaly 

By virtue of the analyticity of the vacuum polarization correlator, the contribution of the 
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hadronic vacuum polarization to a M can be calculated via the dispersion integral |2T 

1 



had 

" ~ 4vr 3 

Ami 



< = 1 — ? / ^had(s)^) • (6) 



Here o"h a d(s) is the total e + e — >hadrons cross section as a function of the cm. energy- 
squared s, and K(s) denotes the QED kernel 

K(s) = x 2 (2 - — ) + (1 + xf ( 1 + \) { ln(l + x) - x + — J + |^4x 2 lnx (7) 



with a; = (1 — + and /3 = (1 — 4m 2 /s) 1 / 2 (see also remarks concerning the 

numerical stability of K(s) in Ref. 0). The function K(s) decreases monotonically with 
increasing s. It gives a strong weight to the low energy part of the integral (||). About 
91% of the total contribution to ajj ad is accumulated at cm. energies y/s below 2.1 GeV 
while 72% of aj^ ad is covered by the two-pion final state which is dominated by the p(770) 
resonance. Recent r data published by the ALEPH Collaboration provide a very precise 
spectrum of the two-pion final state as well as new input for the more controversial four- 
pion final states. This new information can significantly improve the a^ ad determination. 

Running of the QED Fine Structure Constant 

In the same spirit we evaluate the hadronic contribution Aa(s) to the renormalized vac- 
uum polarization function II' (s) which governs the running of the electromagnetic fine 
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structure constant a(s). For the spin 1 photon, II' (s) is given by the Fourier trans- 
form of the contraction of the electromagnetic currents j£ m (s) in the vacuum (q^q u — 
q 2 g^) n;(g 2 ) = i J d*x e iqx (0\T(jg m (x)j» m (0))\0). With Aa(s)=-4vra Re [W^s) - W 7 (0) 
one has 

, ^ «(°) 
a(s) — 



1 - Afl(s) ' 

where 47ra(0) is the square of the electron charge in the long- wavelength Thomson limit. 
The contribution Aa(s) can naturally be subdivided in a leptonic and a hadronic part. 
Furthermore, at s — M| it is appropriate to separate the leading vacuum polarization 
contribution involving the five light quarks u, d, s, c, b from the top quark contribution 
since the latter cannot be calculated in the light fermion approximation. 
The leading order leptonic contribution is given by 

W^) = ^E(^- I) =314.2X10-. (9) 

Using analyticity and unitarity, the dispersion integral for the contribution from the light 
quark hadronic vacuum polarization Aa^j (M|) reads [22 



Acfi(M z 2 ) = -—i-Re I ds—^p— , (10) 
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where a(s) = 167r 2 a 2 (s)/s ■ Imll' (s) from the optical theorem, and Imll^ stands for the 
absorptive part of the hadronic vacuum polarization correlator. In contrast to a^ ad , the 
integration kernel favours cross sections at higher masses. Hence, the improvement when 
including r data is expected to be small. 

The top quark contribution can be calculated using the next-to-next-to-leading order 
prediction of the total inclusive cross section ratio R, defined as 

(Ttot(e + e~ -> hadrons) 3s 



from perturbative QCD |23|, || : 
R 



L pert I 




H h ri — +r 2 — 

7T \ 7T / \ 7T 



(12) 



where n = 1.9857 - 0.1153 n/, r 2 = -6.6368 - 1.2001 n/ - 1.2395(E/ <5/) 2 /3 £/ and 
n/ is the number of involved quark flavours. The evaluation of the integral fllPf) with 
^to P = 175 GeV and the running strong coupling constant fixed at a s (M|) = 0.121 
yields Aa top (M|) = -0.6 x 10" 4 . 

Using Aaf, 5 a ) d (M|)= (280 ± 7) x 10" 4 §, one obtains 

a -1 (M|) = 128.902 ± 0.090 . (13) 

Again, the error is dominated by the hadronic vacuum polarization part that is not 
calculable within perturbative QCD. 



5 



2 The Integration Procedure 



The information used for the evaluation of the integrals (19) and (|H]) comes mainly from di- 
rect measurements of the cross sections in e + e~ annihilation and via CVC from r spectral 
functions. The integrals themselves are evaluated using the trapezoidal rule, i.e., combin- 
ing adjacent measurement points by straight lines. Even if this method is straightforward 
and free from theoretical constraints (other than CVC in the r case) , its numerical calcu- 
lation requires special care. The combination of measurements from different experiments 
taking into account correlations - both inside each data set and between different exper- 
iments are subjected to additional studies and discussions. 



Averaging Data from different Experiments 

In order to exploit the maximum information, we combine weighted measurements of 
different experiments at a given energy instead of calculating separately the integrals for 
every experiment and finally averaging themF]. Generally, if different measurements at a 
given cm. energy show inconsistencies, i.e., their x 2 per degree of freedom (dof) is larger 



than one, we rescale the error of their weighted average with yx 2 /dof. 

The solution of the averaging problem is found using a correlated \ 2 minimization. 
We define 

A'exp N n 

x 2 = E E (*? - (q;-)" 1 (4 - ki) , (14) 

n=l i,j=l 

where is the ith cross section measurement of the nth experiment in a given final state, 
C™j is the covariance between the ith and the jth measurement and k{ is the unknown 
distribution. The covariance matrix C n is defined as 



f (A£ stat ) 2 + (A« sys ) 2 for i 
' ' ~ \ A£ sys ■ A£ sys for i j 



i,j = l,...,N n , (15) 



where A" stat (A" sys ) denotes the statistical (systematic) error of x™. The systematic errors 
of the e + e~ annihilation measurements are essentially due to luminosity and efficiency 
uncertainties. The minimum condition dx 2 /dki = 0, Vz leads to the linear equation 
problem 

^exp N„ 

EE^i-y^r 1 = °> i = h...,N n . (i6) 

n=l j=l 



4 One could imagine two experiments a and b, each with two independent measurements ai(£'i)± Aai, 
02 (£2) ± Aa 2 and b x {E x ) ± A&i, b 2 {E 2 ) ± Ab 2 at energies E x ^ E 2 . Setting Aai = Aa 2 = A&i = Ab 2 
leads to identical errors in both integration methods. However, non-symmetric errors as, e.g., Aai/2 = 
2Aa 2 = 2A6i = Ab 2 /2 propagate a 53% larger uncertainty when calculating independently the sum over 
the points (trivial integration) of the experiments a, b and averaging afterwards rather than averaging 
(cti, bi) and (a 2 , b 2 ) first, i.e., keeping the energy information of the respective points in the average. 



6 



The inverse covariance between the solutions ki, kj is the sum over the inverse 
covariances of every experiment 

c- x = E^r 1 - ( 17 ) 

n=l 

Correlations between Experiments 

Eq. ([IT]) provides the covariance matrix needed for the error propagation when calculating 
the integrals over the solutions fc, from Eq. (|16j). Up to this point, only contains corre- 
lations between the systematic uncertainties within the same experiment. However, due 
to commonly used simulation techniques for acceptance and luminosity determinations as 
well as state-of-the-art calculations of radiative corrections, systematic correlations from 
one experiment to another cannot be excluded. It is clearly a difficult task to reason- 
ably estimate the amount of such correlations as they depend on the capabilities of the 
contributing experiments and one's theoretical understanding of the dynamics of the re- 
spective final states. In general, one can state that in older experiments, where only parts 
of the total solid angle were covered by the detector acceptance, individual experimental 
limitations should dominate the systematic uncertainties. Potentially common system- 
atics, such as radiative corrections or efficiency, acceptance and luminosity calculations 
based on the Monte Carlo simulation, play only minor roles. The correlations between 
systematic errors below 2 GeV cm. energy are therefore estimated to be between 10% 
and 30%, with the exception of the 7r + 7r~ final state, where we impose a 40% correlation 
due to the easier experimental situation and the better knowledge of the dynamics which 
leads to non-negligible systematic contributions from the uncertainties of the radiative 
corrections. At energies above 2.1 GeV the experiments measured the total inclusive 
cross section ratio R. Between 2.1 and 3.1 GeV, individual technical problems dominate 
the systematic uncertainties. At higher energies, new experiments provide nearly full ge- 
ometrical acceptance which decreases the uncertainty of efficiency estimations. Radiative 
corrections as well as theoretical errors of the luminosity determination give important 
contributions to the final systematic errors quoted by the experiments. We therefore es- 
timate the correlations between the systematic errors of the experiments to be negligible 
between 2 GeV and 3 GeV, 20% between 3 GeV and 10 GeV, and 30% above 10 GeV. 

These correlation coefficients are added to all those entries of from Eq. (|i"7|) which 
involve two different experiments. 

Inclusion of r Data 

In this analysis we include additional data from r decays into two- and four-pion final 
states 5 : r~ ^ tt^tt u T , t~ ^ tt^Stt v t and r~ — > 2tt~7t + 7t u t , recently published by the 
ALEPH Collaboration Assuming isospin invariance to hold, the corresponding e + e~ 
isovector cross sections are calculated via the CVC relations 

i=i 47T « 2 / 1B N 

S 

5 Throughout this paper, charge conjugate states are implied. 
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o- e+e -_ Mr+7r - 7r+7r - = 2 ■ —v ltW - 3n o UT , (19) 



2 



i=i 4:7ia 2 r . . . 

CX e+e -_ 7r+7r - 7r 7r = — ^— K^Tr-Tr+Tr ^ ~ ^1, 7T" 37r<V J ■ ( 20 ) 

The r spectral functions v\ are given as binned continous distributions of the mass- 
squared s. In each bin i, the spectral function v\(si) contains the (normalized) invariant 
mass spectrum ANi/N, integrated over the bin width dsi. On the contrary e + e~ cross 
sections are measured at discrete energy settings. To each e + e~ measurement is associated 
a t cross section value obtained by interpolating between adjacent bins. This interpo- 
lation takes into account the correlations between the bins and is achieved imposing a 
functional form obtained from a fit to Breit-Wigner resonances |§ using the Gounaris- 
Sakurai parametrization [25|| . However, the fit function is renormalized in each bin so that 
its integral over the width of each bin corresponds to the measured CVC cross section for 
that bin. All the data points - r, e + e~ and interpolated r values - are injected with their 
correlations into Eqs. (pl|) and ([H 



Due to the high bin-to-bin correlations of the r data and the significant normalization 
uncertainty coming from the r hadronic branching ratios, biases of the least-square mini- 
mization cannot be excluded. We therefore calculate the average solution twice, i.e., 



with and without correlations, use the mean value of both integrations as the result and 
add half of the total difference as additional systematic error. This is done in all cases 
where r data are involved. The effect amounts to about 10% of the total error. 



Evaluation of the Integral 

The procedure described above provides the weighted average and the covariance of the 
cross sections from different experiments contributing to a certain final state in a given 
range of cm. energies. We now apply the trapezoidal rule. In order to perform the 
integrations (^j and (|1U|), we subdivide the integration range in sufficiently small energy 
steps and calculate for each of these steps the corresponding covariance (where additional 
correlations induced by the trapezoidal rule has to be taken into account). This procedure 
yields error envelopes between adjacent measurements as depicted by the shaded bands 
in Figs. H and ||. 



3 Radiative Corrections 

Higher order radiative corrections bias the measurements of both cross sections in e + e~ 
annihilation and invariant mass spectra from r hadronic decays. The e + e~ experiments 
generally correct the measured cross section for QED effects, including bremsstrahlung, 
vacuum polarization and higher order self-energy graphs (see references in ) . Following 
the prescription of Ref . , we have multiplied all inclusive cross section measurements R 
at masses below the J/i/j resonance by the (small) correction factor (l+Aai ep (s))(a/a(s)) 2 
in order to account for the missing correction for hadronic vacuum polarization. 

In r decays, final state radiation from the r itself or from its decay products can 
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influence the invariant mass measurement. Due to the high mass of the r lepton, the 
bremsstrahlung graph is largely suppressed. Both types of radiation are included in the 
Monte Carlo simulation program KORALZ p6f , used by ALEPH to unfold the measured 



distributions from detector resolution and physical (higher order) effects. Even if the 
actual frequency with which final state radiation occurs or if its energy was not well 
simulated in the Monte Carlo, reconstructed photons found to originate from radiation of 
the r decay |27| are included in the invariant mass determination, and thus do not bias 



the measurement. 

Electroweak radiative corrections are applied through the CVC correction factor 5*ew 
defined in Eq. (El). 



4 The Origin of the Data 

The exclusive low energy e + e~ cross sections have mainly been measured by experiments 
working at e + e~ colliders in Novosibirsk and Orsay. Due to the high hadron multiplicity 

at energies above 2.5 3.1 GeV, the exclusive measurement of the respective hadronic 

final states is not practicable. Consequently, the experiments at the high energy colliders 
DORIS and PETRA (DESY) and PEP (SLAC) have measured the total inclusive cross 
section ratio R. 

We give in the following a compilation of the data used in this analysis: 



The e + e 



+ " r measurements are taken from OLYA 



cmd m 



-> 7T 1 7T 

DM1 and DM2 |3 
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3, TOF g, NA7 |1| 
In addition, we use r data from ALEPH H nor 



malized to the world average branching ratio B(r — > n 7r° v T ) = (25.24 ± 0.16)% [1 1 
According to Eq. (|i~8l) , r data provide only the dominant isovector part of the to 



tal two-pion cross section. A correction due to the small isospin-violating isoscalar 
uj — > 7r + 7r~ final state, which interferes with the isovector amplitude, is applied. A 
small correction for the missing, i.e., unmeasured decay modes p — > 7r°7 (only for 
e + e~ data) and p — > 7/7, is added. 

- The reaction e + e~^ 7r + 7r~7r° is dominated by uj and <fi intermediate resonances. In 
the peak region of these resonances we use analytic parametrizations of the cross 
sections. The non-resonant data are taken from ND [3jj], M2N [|35|] , M3N |36|, 
DM1 WR and DM2 |38||. Corrections for the missing uj and <fi decay modes are 



applied. 

The e + e 
DM2 M 



7T n TC 



-7r°7r° 



41 



data are available from OLYA 
According to Eq. 



ND 151, M2N 1351 



and M3N 36 



|20f) , a linear combination of 
both four-pion r decay channels measured by ALEPH H] connects the corresponding 
spectral functions with the above e + e~ final state. We use the two branching ratios 

(4.25 ±0.09)% and B(r" -> Tr^Tr v T ) = (1.14 ± 0.14)% P 



B(t~ ->27r-7r 
as an appropriate normalization. 



The reaction e + e — > ujtt is mainly reconstructed in the 7t + tt 7r°7r° final state. It 
was studied by the collaborations ND [34] and DM2 [42]. Corrections for the missing 
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u decay modes are applied. 



The e + e — > tt + tt tt + tt final state was studied by the experiments OLYA [3S 
ND H, MEA g, CMD ||, DM1 ||, 0, DM2 
The corresponding spectral function from r" 
provided by ALEPH §. 



J |42], |4J and M3N 

7r _ 37r°z/ T (according to Eq. (|T9|)) is 



The e + e~^ 7r + 7r-7r + 7r-7r final state is taken from M3N g§] and CMD g§. The 
other five-pion mode e + e~ — > 7r + 7r~37r° can be accounted for using the rigourous 
isospin relation a^+^-^o = 0.5 x cr 7r+7r - 7r+7r - 7r o. 



For the reaction e + e _ — > ton + n~ , measured by the groups DM1 [[4(J an d DM2 [[45 
a correction for u decays other than into three pions which appear in the five-pion 
final state is applied. 

The e + e~ — > tt + 7t~7] data were studied by ND |4f and DM2 |g8). We subtract from 



the cross section the contributions which are already counted in the 7r H 
and 7r + 7T~37r° finale states. 



7T 7T + 7T 7T 







The cross sections of the six-pion final states 37r + 37r and 2ir + 2ir 27r° were measured 
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by DM1 [49|, M3N [36] CMD 



41]. In Ref. M an upper limit for 



and DM2 

the unknown n~7i + 4ti cross section of o" 7r + 7r _ 47r o < (3/2) x cr 2vr - 27r + 27r o — (9/24) x 
o"3 7r - 37r + was derived using isospin constraints. We take half of this upper limit as 
the estimated contribution, with an error of 100%. 



The e+e^^K+K" and e+e' 



K° S K° 



cross sections are taken from OLYA Q 



DM1 H and DM2 [|2 



The reactions e + e~^KgK + 7r - and e + e~^K + K~7r were studied by DM1 p3" , |54 | 
and DM2 [40]. Using isospin symmetry the cross section of the final state KgR^vr 
is obtained from the relation OkOko^o = uk+k-tt - 



The inclusive reaction e + e — >Kg+X was analyzed by DM1 [ 55| . After subtract- 
ing from its cross section the separately measured contributions of the final states 



KgK+TT 



K° S K- 



7T 



and KqK?7r°, it still includes the modes KcK 



7T 







and KgK" 



7T 



r0 





s- lv s 



7T 7T 



7T 7T 



0X^0/ 



7r". With the assumption that the cross sections for the 
processes e T e" — >.K.".K."(7r7r / ) and e + e~ ^K + K~(7T7r)° are equal, one can summarize 
the total KKtttt contribution as twice the above corrected K s +A cross section. A 
reasonable estimate of the systematic uncertainty, implied by the assumption made, 
is obtained by taking the cross section for the channel K + K~7r + 7r~ measured by 
DM1 E3I and DM2 ffl. 



At higher energy the total cross section ratio R is measured inclusively. We use the 

MARK I [§7 
LENA !>Ii, Crystal Ball HB3, MD-1 

m 



data provided by the experiments 772 
PLUTO 0, LENA frf], Crystal Ball §2 
MARK- J §§, TASSO 0, CLEO @ 



§, DASP m 



m 



DELCO ^ 

CELLO @, JADE |S 
and MAC 0. Above 3.5 GeV 



[63 



CUSB 

the measurements of the MARK I Collaboration are significantly higher than those 
from LENA, PLUTO and Crystal Ball. In addition, the QCD prediction of R, which 
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should be reliable in this energy regime, favours lower values. In agreement with 
Ref. 0, we neglect MARK I data above this energy threshold. 

Although small, the enhancement of the cross section due to 7-Z interference is 
corrected for cm. energies above the J ftp mass. We use a factorial ansatz according 
to Ref. [[Ti], ||, yielding a negligible contribution to a^ ad and a —0.30 x 10~ 4 shift of 



5 Analytical Contributions 

In some energy regions where data information is scarce and/or reliable theoretical pre- 
dictions are available, we use analytical contributions to extend the experimental integral. 



The 7r + 7r Threshold Region 

To overcome the lack of data at threshold energies, a second order expansion obtained from 
Chiral Perturbation Theory J72| is used as a description of the pion form factor F n (which 
is connected with the two-pion cross section via the expression {F^ 2 = 3sa W7T /(TTa 2 (l — 
Aml/sfl 2 ): 

F. ChPT * l + Ur\s + c w s 2 + 0(s') . (21) 


Exploiting precise results from space-like data [ff^l , the pion charge radius-squared (r 2 ) n = 
(0.431 ± 0.026) fm 2 and the coefficient = (3.2 ± 1.0) GeV" 4 from Eq. (gTJ) have recently 
been determined WM by means of a simultaneous fit. 



Narrow Resonances 

The e + e~ annihilation cross section involves narrow resonances such as the a;(782) and 
0(1020) at low energies, the J/ip and T resonances at the cc and bb quark thresholds, 
respectively, as well as their excited spectroscopic states. It is safe to parametrize these 
states using relativistic Breit-Wigner resonance shapes with an s-dependent width. We 
use the formulae given in Ref. || . The physical input values of the parametrizations and 



their errors are taken from Ref. [|TTJ . The total parametrization errors are then calculated 
by gaussian error propagation. 



High Energy Tail 

At energies sufficiently above the T resonance family, the perturbative QCD prediction 
of R with five active quarks is supposed to be reliable. In agreement with Ref. we use 
R peTt (s) from Eq. (P) for yfi > 40 GeV. 



6 Results 

We evaluate the integrals @ and ([U]) exclusively, i.e., for every contributing final state, 
up to the cm. energy of 2.125 GeV. Even if some particular modes have been measured 
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up to somewhat higher energies, we have to worry about unmeasured exclusive modes 
and therefore use the total R measurement above this threshold. Both energy regions 
are assumed to be uncorrelated. Because the contributions of the exclusive channels 
at low energy are simply summed up, we have to estimate their respective covariances 
when propagating the error: in general, unmeasured final states whose contributions 
are deduced from measured ones via isospin are set to be 100% correlated with these. 
Also different detectors performing the same measurement are correlated through the 
sharing of commonly used simulation techniques to calculate acceptance and selection 
efficiency which depend on the assumed underlying physical dynamics. Contributions 
from resonances that are analytic are globally assumed to have 20% correlations due to 
modelling uncertainties. Between purely measured final states we have estimated the 
correlations depending on the number of common experiments that contribute to their 
measurements and on the common energy region, as well as according to the relative 
importance of their statistical and systematic errors. In general our estimation yields a 
correlation between 10% and 20%. This treatment is different from that of Ref. J75|| where 
a 100% correlation was assumed. 

As described in Section p], corrections to the charged p ± width have to be applied to 
account for small CVC-violating effects. The magnitude of the width difference Eq. (^) 
translated into a^ ad and Aa[^ d (M|) is evaluated using a parametrization of the p line 
shape based on vector resonances ||. We obtain the additive corrections 

Sa had = _( L3±2 .0) x 10~ 10 
SAa { ^ d (Mi) = -(0.09 ± 0.12) x 10" 4 (22) 

for the r~ — > 7t~7t° v T spectral function which is applied in the present analysis. Corrections 
from the higher mass resonances p(1450), p(1700) are expected to be negligible. 

The two- and four-pion cross sections (incl. the r contribution) in different energy 
regions are depicted in Figs. |2| and |3|. The bands are the results within (diagonal) error- 
envelopes of the averaging procedure and the application of the trapezoidal rule described 
in Section |2|. 



Lowest Order Hadronic Contributions 



The results of the exclusive contributions to a^ ad and Aa|^ d (M|) are presented in Ta 



ble 6. After the inclusion of the r data, the error of aj^ ad is dominated much less by the 
uncertainty of the two-pion contribution. 

Other important sources are the contradictory 7r + 7r~27r° data, as well as the unmea- 
sured 7T + 7r - 47r° final state. In the latter case, limits can be set only by using very 
conservative isospin arguments. As shown in Ref. |J by means of a decomposition in 
orthogonal isospin classes (Pais [76[]) the large upper limit comes from the assumption of 



a dominant a^n class accompanied by a vanishing 0321 contribution. Both classes occur 
in 7r + 7r~7r + 7r~27T°, while none contributes to 7r + 7r~7r + 7r~7r + 7r~ and only oaa_\ is part of 
7T + 7r~47r . The measured cross section of 7r + 7r~7r + 7r~27r° is clearly higher than the cor- 
responding 7T + 7r~7r + 7r~7r + 7r~ final state hence guaranteeing a leading contribution from 
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Figure 2: Two-pion cross section as a function of the cm. energy- squared. The band 
represents the result of the averaging procedure described in the text within its diagonal 
errors. The lower left hand plot shows the chiral expansion of the two-pion cross section 
obtained from expression (\2l\). 
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Final states 


ajf (xlO iU ) 


Mad( M l) ( Xl ° ) 


Energy (GeV) 


7r + 7r _ threshold 


2.30 ± 0.05 


0.04 ± 0.00 


; rj 

Am{- 0-320 


7T + 7T~ 


495.86 ± 12.46 


34.01 ± 0.87 


0.320-2.125 


7r + 7r~ (incl. r data) 


500.81 ± 6.03 


34.31 ± 0.38 


0.320-2.125 


p(7r 7 + ?7 7) W 


0.30 ± 0.05 


0.02 ± 0.01 


0.298-2.125 


a; 


37.09 ± 1.07 


2.97 ± 0.09 


0.420-0.810 


a; — ► 7T7, neutrals 1 > 


0.03 ± 0.01 


< 0.01 


0.810-2.125 




39.23 ± 0.94 


5.18 ± 0.12 


1.000-1.055 


$ — > 777, 7r u 7 ^ 


0.09 ± 0.01 


0.01 ± 0.00 


1.055-2.125 


7r + 7r~7r u (below <£) 


4.12 ± 0.41 


0.42 ± 0.04 


0.810-1.000 


71^ 7r 7r (above <E>) 


1.90 ± 0.72 


0.46 ± 0.26 


1.055-2.125 


7T+7T 27T U 


21.41 ± 2.36 


5.82 ± 0.63 


0.910-2.125 


7r + 7T _ 27r° (incl. r data) 


22.26 ± 1.53 


6.16 ± 0.49 


0.897-2.125 


W7r°(u; — > 7T7, neutr.) W 


0.88 ± 0.11 


0.18 ± 0.02 


0.930-2.125 




15.90 ± 1.34 


4.61 ± 0.39 


0.983-2.125 


7r + 7r _ 7r + 7r _ (incl. r data) 


16.50 ± 0.98 


4.76 ± 0.31 


0.794-2.125 


1 | A 


4.02 ± 0.51 


1.51 ± 0.20 


1.019-2.125 


7T+7T 37T U ^ 


2.01 ± 0.26 


0.75 ± 0.10 


1.019-2.125 


ujtt^tt [co — ► 7T7, neutr.) 1 ; 


0.07 ± 0.02 


0.03 ± 0.01 


1.340-2.125 


7T + 7T ~ 7T + 7T ™ 7T + 7T ~~ 


0.47 ± 0.14 


0.19 ± 0.04 


1.350-2.125 


7T 7T 7T^7T 27T 


3.32 ± 0.36 


1.35 ± 0.14 


1.350-2.125 


7T+7T 47T U ^ 


2.40 ± 2.40 


0.98 ± 0.98 


1.350-2.125 




0.51 ± 0.14 


0.16 ± 0.05 


1.075-2.125 


K+K" 


4.30 ± 0.58 


0.85 ± 0.10 


1.055-2.055 




1.20 ± 0.42 


0.23 ± 0.08 


1.090-2.125 


KgK + 7r (+ K^K -k +{Z >) 


2.04 ± 0.36 


0.70 ± 0.12 


1.340-2.125 




0.42 ± 0.29 


Air i r\ i n 

0.15 ± 0.10 


1.440-2.125 


t^Ot^O^O (2) 


n /I o _i_ n on 
U.4z ± U.zy 


n i c ! n in 
U.lO ± U.1U 


1 A AC\ T IOC 

1.44U — Z.lzo 


KK7T7T (all modes) 


4.52 ± 1.65 


1.82 ± 0.66 


1.441-2.125 


J/V>(1S,2S,3770) 


8.04 ± 0.52 


9.97 ± 0.68 


3.096-3.800 


T(1S,2S,3S,4S,10860,11020) 


0.10 ± 0.01 


1.18 ± 0.08 


9.460-11.20 


R 


41.64 ± 3.61 


164.31 ± 5.59 


2.125-40.0 


R (perturbative) W 


0.16 ± 0.00 


42.82 ± 0.10 


40.0-oo 


J2 (e + e~ — >hadrons) 


695.0 ± 15.0 


280.9 ± 6.3 


4m^ -oo 


(e + e~ — >hadrons) (incl. r data) 


701.1 ± 9.4 


281.7 ± 6.2 


4m^ -oo 



1 Correction for missing modes (see text). 

2 Deduced from isospin relations (see text). 

Without contribution from r\ — > ir + ir~ir a and r] — ► 3tt° . 
4 Values are taken from 0j . 

Table 2: Summary of the ajj ad and Aa^ d (M|) contributions from e + e~ annihilation and 
t decays. The line threshold" contains the results from the integral over expres- 

sion $n\ ) at threshold energies. 
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Figure 3: Four-pion cross section as a function of the cm. energy-squared. The band 
represents the result of the averaging procedure described in the text within its diagonal 
errors. The right-hand plot shows additionally the e + e~— > ujtt amplitude (small points). 



one of the classes mentioned if isospin invariance holds. Since those classes correspond 
to eigenstates, a resonance analysis of the measured six-pion data would reveal impor- 
tant properties of the class structure of the respective modes which thus could give more 
constraining isospin bounds. 

Another large uncertainty comes from the KKirn final states. The measurement of 
the K + K~7r + 7r~ mode alone does not allow one to calculate isospin bounds for all possible 
contributions. Fortunately, it is possible to extract the complete KKtctt contribution on 
the basis of a DM1 measurement of the inclusive channel K s +X [55|| . Nevertheless, large 
experimental uncertainties prevent a precise determination of the corresponding integrals. 

The last important error source, especially for Aa^ d (M|), comes from the integral 
over the measured high energy inclusive cross section ratio R. The reliability of the QCD 
perturbative expansion for energies sufficiently above the still unpredictable resonance 
phenomena has been proven in many cases (see, e.g., a s measurements from different 
energy scales at LEP, HERA and from r decays). Thus, theoretical input at energies 
lower than 40 GeV should be reliable and could significantly help to reduce the integration 
uncertainties [[79|, [80] . However, this has not been used in the present analysis which relies 
on experimental data as far as possible. 

The squared contributions of the various final states and energy regimes to the errors 
of aj^ ad and Aa(M|) are depicted in Fig. §. Only the results after the inclusion of r data 
are shown. 
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A 



Figure 4: Quadratic contribution of the various error sources to aj^ ad (left hand plot) and 

Aaj^ d (Mf) (right hand plot) after the inclusion of t data. In the energy region 0-2.1 GeV 
we include all exclusive contributions that are not given separately. 



We obtain for the lowest order hadronic vacuum polarization diagram of the muonic 
anomalous magnetic moment the contributions 



a 



had _ (ac\K n i icn\ v 1 n-10 
had _ fi-rn-l 1 _l_ n A\ w 1 n-10 



(695.0 ± 15.0) x 10~ 10 (e + e~ data only) 



fl had = (yon ± 9 4) x iq-lv (combined e + e" and r data) 

and for the running of a at M| 

Aa h J d (M|) = (280.9 ± 6.3) x 10" 4 (e + e~ data only) 
AaSi(Aff) = (281.7 ± 6.2) x 10" 4 (combined e + e" and r data) . 

Fig. |5| shows a compilation of published results. The inclusion of the new r data 
yields a large improvement in the precision of the aj^ ad determination. The difference 
in aj^ ad between the exclusive e + e~ analyses of and this work is mainly due to a 
disagreement in the two-pion integral where we obtain significantly lower values. In 
addition, differences in the handling of unmeasured modes generate inconsistent results. 
The results of Ref. cannot easily be compared to the newer ones as the data set did 
not include the recent DM2 results. The differences in the final errors of a^ ad in the 



exclusive e + e~ analysis of this work compared to Ref. [y, 75 1 is mainly caused by different 
techniques in the handling of the data and their errors. The detailed study of the origin 
of correlations, their propagation, as well as the rigourous use of isospin constraints to 
bound unmeasured modes yield slightly smaller errors here. As expected, the gain in the 
precision of aj^ ad coming from r data is significant whereas it is very small for Aa^M^) 
since the dominant contributions and uncertainties come from energies above the r mass. 
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Kinoshita et al. [7] 
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Figure 5: Comparison of estimates of a^ ad (lowest) and Aa|^ d (M| 



Higher Order Contributions 

In the famous paper of Kinoshita et al. 0, higher order contributions to the muon hadronic 
vacuum polarization graph, such as additional lepton or quark loops inserted in the dia- 
gram of Fig. [I] and the so-called light-by- light scattering graph, have been evaluated. The 
latter has been re-computed by different groups obtaining (—5.2 ± 1.8) x 10~ 10 pl| and 
(—9.2 ± 3.2) x 10~ 10 |33] with large uncertainties compared to the designed experimental 
accuracy of Aa /t ~ 4.0 x 10~ 10 of the forthcoming BNL experiment. We use the value of 
(—6.2 ± 4.0) x 10~ 10 in the following with an enlarged error to account for inconsistencies. 

The calculation of the higher order 0(a/7r) 3 loop diagrams is accomplished and partly 
corrected in a recent work where second order kernel functions K^(s) are provided. 
These are used to calculate the corresponding contributions in the same spirit as the dom- 
inant lowest order graph by virtue of the dispersion integral (|B|). The numerical evaluation 
in Ref. |84| was performed on the basis of the data sample used by Ref. ||. We repeat 
this exercise here in order to check the consistency of the results. For the contribution of 
diagrams with additional photon exchanges, e.g., the fourth order muon vertex correction, 
we use the kernel labeled K {2a \s) in fg§ and obtain a% a) = (-20.9 ±0.4) x 10" 10 . The 
diagrams with an electron loop inserted in one of the photon lines of Fig. [j] (kernel K ( 2b ' (s) 
in Ref. |83]]) contribute to = (10.6 ±0.2) x 10~ 10 , where the asymptotic expansion, 



the analytical and numerical solutions provided in Ref. |84| lead to very similar results. 
Finally, the insertion of two hadronic loops in the muon vertex correction graph (kernel 
K( 2c \s) in Ref. |3J) results in af^ = (0.27 ±0.01) x 10~ 10 . The contributions a {2a > b > c) are 
found to be in agreement with Ref. |84|| . All higher order results given here are computed 
from the e + e~ data set only. 

The compilation of the hadronic higher order parts (including light-by-light scattering) 



yields < ad [(a/7r 



■16.2 ± 4.0) x 10 



-10 
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Results for and cx(M^) 

Collecting all contributions, we obtain for the anomalous magnetic moment of the muon 

a M = (11659164.5 ± 15.6) x KT 10 (e + e" data only) 

a M = (11659170.6 ± 10.2) x 10~ 10 (combined e + e~ and r data) , 

where the errors of the lowest order calculation aj^ ad and 

a (2a,b,c) are 

added linearly. 

The inverse of the fine structure constant at M| is found to be 

a' 1 (Ml) = 128.889 ± 0.087 (e + e" data only) 

a" 1 [Ml) = 128.878 ± 0.085 (combined e + e" and r data) . 

One may use the latter (combined) result for a(M|) to improve the constraint on the 
mass of the standard model Higgs boson Mmggs inferred from a global electroweak fit. 



This is done by utilizing current available electroweak data [85], |86| and the ZFITTER 
electroweak library [j^Tf] . It requires M z , m top , and a(M|) as input parameters, which are 
allowed to vary within their experimental accuracies. The additional parameters MHi ggs 
and the strong coupling constant at Mf , a s (M|), are freely adjusted in the fit. We obtain 
a s (Ml) = 0.1201 ± 0.0033 which is in perfect agreement with the experimental value 



of 0.122 ± 0.006 |SB| from the analyses of QCD observables in hadronic Z decays at LEP. 
The fitted Higgs boson mass is 138i^ 7 GeV, compared to 149t^ 8 GeV when usine; the 
previous value of a(M|) from Eq. (|13|) . An additional error of 50 GeV should be added 
to account for theoretical uncertainties f87j| . 

Fig. P depicts the variation of x 2 as a function of the Higgs boson mass for the new 
and previously used values of «(M|) (the latter taken from ||). We obtain an upper limit 
for M ffig gs of 516 GeV at 95% CL. 



7 Conclusions 

We have reevaluated the hadronic vacuum polarization contribution to (g — 2) of the muon 
and to the running of the QED fine structure constant a(s) at s = M§. We used new data 
from r decays, recently published by the ALEPH Collaboration, in addition to slightly 
enlarged e + e~ annihilation cross section data sets in order to improve the precision of 
the corresponding integrals. Our results are, to lowest order, a^ ad = (701.1 ±9.3) xlO -10 

yielding a M = (11659 170.6 ± 10.2) x lO" 10 and Aag^Mf) = (281.7 ± 6.2) xlO" 4 , propa- 
gating a^ 1 (0) to a _1 (M|) = (128.878 ± 0.085). The improvement coming from r data is 
small in the Aa£ } d (M|) case which is dominated by high energy contributions. However, 
it causes a 37% reduction in the error on a^ ad . 

In the near future, new low energy e + e~ annihilation data are expected to be pro- 
duced by the CMD-2 Collaboration at Novosibirsk. In addition, new results for the 
r~ ^7t~7t°v t spectral function with a precision comparable to the ALEPH data were 
recently presented by the CLEO Collaboration |89"1 |. Significant improvement is also ex- 



pected from energy scans at the future high-luminosity e + e collider DA$NE in Frascati. 



18 




M m ggs (GeV) 

Figure 6: Constraint fit results for the previous and the new value o/a(M|) as a function 
of the Higgs mass. 
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